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ABSTRACT 


A  theorem  is  proved  to  the  effect  that  a  wave  function  for 
a  set  of  N  spins,  which  is  a  product  of  single- spin  wave 
functions  and  which  is  an  eigenstate  of  the  square  of  the  total 
spin  S,  must  be  a  state  with  the  maximum  possible  value  of  S 
and  of  S  ,  z  being  an  arbitrary  direction.  This  theorem  has 
been  applied  in  a  separate  work  by  the  authors  to  show  strik¬ 
ingly  that  the  imposition  of  symmetry  restrictions  of  a  common 
type  on  an  approximate  wave  function  can  lead  to  a  very  poor 
physical  description. 
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A  THEOREM  ON  SPIN-EIGENFUNCTIONS 


I  INTRODUCTION 

In  a  paper  by  the  authors  on  Hartree-Fock  (HF)  theory, 

[ Phys .  Rev.  156,  1  (1967)],  an  analogy  is  drawn  between  the  HF 

approximation  to  a  ground  eigenfunction  of  a  Hamiltonian  H 

representing  a  system  of  electrons,  and  the  Hartree  approximation 

to  the  ground  eigenstate  of  a  Hamiltonian  H  representing  a 

s 

system  of  spins .  An  example  is  then  presented  which  shows  strik¬ 
ingly  that  making  symmetry-restrictions  on  one's  already  other¬ 
wise  restricted  wave  function  (a  rather  common  procedure)  can 
lead  to  an  extremely  poor  description. 

Namely,  let  H  be  an  isotropic  Heisenberg  Hamiltonian 
s 

-Z  j  X  with  exchange  parameters  ^  chosen  in  such  a  way 
that  the  Hamiltonian  represents  an  ant if erromagnet .  The  are 
the  individual  spin  operators  satisfying  by  definition 

— ►  2 

=  S  (S+l) ,  independent  of  i.  The  Hamiltonian  commutes  with 
—►2  — ►  N  — ► 

S  ,  where  S  =  2  s.  .  The  exact  energy  eigenstates  can,  there- 

1—  _L  1 

_+2 

fore,  be  chosen  to  be  eigenstates  of  S  also,  and  it  is  often 
helpful  in  practice  to  require  this  when  dealing  with  exact 
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eigenstates .  This  suggests  (following  others)  that  we  impose 

the  same  requirement  on  the  Hartree  approximation  to  the  wave 

function.  The  latter,  by  definition,  is  of  product  form 

¥  =  II  .<*> .  ,  where  <J>  .  is  a  single-spin  state.  Let  us  therefore 

2 

require  that  Y  be  an  eigenstate  of  S  .  The  problem  then,  in 

this  restricted  Hartree  theory,  is  to  determine  a  lowest  energy 

—►2 

product  wave  function  ¥  that  is  an  eigenstate  of  S  .  (This  is 
completely  analogous  to  a  type  of  restriction  conventionally 
used  in  symmetry- restricted  HF  theory) . 

We  then  make  use  of  a  theorem ;  a  product  ¥  of  spin-s  spin 

—*■2 

functions  which  is  an  eigenfunction  of  S  has  maximum  multi¬ 
plicity  and  the  component  of  S  in  some  direction  has  the  maximum 

possible  value.  According  to  this  theorem,  the  only  ¥'s  of 

—►2 

product  form  which  are  eigenstates  of  S  are  ferromagnetic 

_*2 

(with  maximum  value  of  S  ).  Thus,  the  symmetry- restricted 
Hartree  approximation  is  the  worst  possible  approximation  in  the 
present  example  of  a  Heisenberg  antiferromagnet . 

In  the  paper  on  Hartree-Fock  theory  mentioned  earlier,  the 
above  theorem  was  stated  without  proof;  the  purpose  of  the 
present  note  is  to  supply  the  required  proof,  which  is  given  in 
the  next  section . 

It  may  be  of  interest  to  the  reader  that  the  theorem  seemed 
to  us  to  be  very  probably  true,  and  we  expected  to  find  a  very 
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simple  proof  without  difficulty.  Unfortunately,  we  have  had  to 
be  satisfied  with  the  following  lengthy  proof. 

II.  STATEMENT  AND  PROOF  OF  THEOREM 

We  first  state  the  theorem  more  precisely  and  then  give  the 
proof . 

Theorem .  Let  ¥  =  0^  (1)02  (2)  .  .  .  .0  (N)  be  a  product  of  spin 

—*•2 

functions  0^  such  that  s^  0^(1)  =  s(s+l)0^(l)  and  <0i|0i>  =  1 
for  i=  1, 2,  . . . ,N .  If 

S-2  1  =  o  (a+l)¥  (1) 

N 

where  S*  =  2  ~s. ,  then  "si  -z  0  .  (1)  =  s0  •  (1) ,  i  =  1,  2,  .  .  .,N  for 

•  i  A  JL  _L  -L 

i=  1 

some  unit  vector  z,  and  hence  a  =  Ns . 

Proof .  Before  giving  the  proof  in  detail  we  sketch  an  out¬ 
line  of  it.  We  first  prove  the  theorem  for  N  =  2,  since  the 
proof  for  general  N  is  relatively  short  and  simple  when  based  on 
the  result  for  N  =  2.  The  proof  for  N  =  2  involves  two  parts 
depending  on  whether  or  not  <02|"s]02>  =  T^e  most  lengthy 

argument  by  far  concerns  the  <02  1*1  02>  =  0  case . 

It  is  convenient  to  define 


Q  = 


2 

i<j 


•  s 


3 


1 

2 


1 

2 


3“  2  -  Ns  (&fl) 


3 


and  the  eigenstate  problem 
Q¥  = 

equivalent  to  (1),  where  X  =  [o(o+l)  -  Ns(s+1)]. 

Consider  (2)  first  for  N  =  2. 

?1*2  01(1)02  (2)  =  X  01(1)02(2) 

Take  the  scalar  product  with  02(2)  to  obtain 

t±  -0*2  01  (1)  =  X  0X  (1) 

where 

^  =  <0il  s‘l0i> 

There  are  two  cases  to  consider:  (a)  f  0  and  (b) 

(a)  ^  0 

From  (4),  0^(1)  is  an  eigenfunction  Xm(l)  of  s*^  *o 2 

V°2  Xm(1)  =  m  Xm(1) 

Expand  02  (2)  in  terms  of  the  Xm 

0O  =  2  b  ,  Xm, 

2  ,  m'  m' 

m' 

and  substitute  into  (3)  using 


(2) 


(3) 


(4) 


(5) 
=  0. 


(6) 


(7) 
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(8) 


Sl'S2  “  SlzS2z  +  2  sl+  s2-  +  2  Sl-  S2+ 


where  s+  =  +  is^.  Obtain ^  using 


s±  Xm  “  g±m  xm±l 


(9) 


1/2 

where  =  [  (s-m)  (s+ith-1) ]  =  9_m_1 


m 


Z,  m'  bm'  Xm(1)V(2)  +  \  9m  Vl(1)  Z,  g-m'  V  V-l 

m'  m' 


(2) 


+  7  g-m  W1'  S,  V  bm-  V+l(2) 

m' 


X  Xm(1)  s ,K'  V  (2) 

m' 


(10) 


Take  the  scalar  product  with  Xm(l)  Xm,  (2),  obtaining  thereby 


(mm' -X)b  .  =  0 
m' 


(id 


This  implies  b^,  =  6m,—  for  some  m,  -s  <  m  <  s,  and  consequently 
X  =  mm.  Substituting  into  (10)  then  yields 


mm 


Xm(1)  Xm(2)  +  2  gm  g-m  Xrtvfl(1)  Xm-1(2) 


+  2  g-m  gH  Xm-l(1)  W2)  =  Xm(D  V2> 


(12) 
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from  which  we  conclude  that 


gm  g  —  =  o 

m  -m 


-m 


g_  =  o 

m 


(13a) 

(13b) 


Equations  (13)  require  either  g  =  g_  =  0  and  hence  m  =  m  =  s 

m  m 


or  g  —  =  g  „ 
-m  -ro 


m  "m 

=  0  and  hence  m  =  m  =  -  s .  Thus,  we  conclude  that 


when  o2  ^  (3)  has  only  solutions  of  the  type 


0  •  =  X 

1  s 


i  =  1,2 


(14) 


where  s,  *z  X  (1)  =  s  X  (1)  with  z  some  unit  vector. 

Is  s 

For  an  alternate  proof  we  note:  once  we  know  <p ^  =  Xmj  then 

we  know  immediately  that  the  only  possible  product  eigenfunctions 

of  S*2of  the  form  X  (1)0~(2)  are  just  x  (1)X_(2)  and 

in  ^  s  s 

*_s(1)X-s(2)  ' 

(b)  a*2  =  0 

In  this  case  X  =  0  (excluding  the  physically  uninteresting 
¥  =  0  case)  and  (3)  reduces  to 


si‘s2  01d)02(2)  =  0 


(15) 


Choosing  an  arbitrary  z-axis  we  expand  0^  and  02 
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(16a) 


0,  =  2  a  Xm 
1  mm 

m 


=  2  b  ,  Xm, 
2  ,  m'  nr 

m' 


(16b) 


where 


s  •  z  Xm  =  m  X 


m 


m 


-s  <  m  <  s 


(17) 


Then,  (15)  becomes,  on  using  (8)  and  (16), 


2  ,  am  bm'  (2)  +  ?  9m  9-m'  am  V  W(1,V-1<2) 

mm'  *>• 


+  -X  g  _  g_,  a  b  X  ,  (1)X_, -i  (2) 

2  -m  m  m  rrr  m- 1  m-fl 


}=° 


(18) 


Taking  the  scalar  product  with  X  (l)Xm, (2)  then  gives 


mm' 


m  m' 


+  -x  gm  n  gm, 

2  m- 1  m 


am-l  ^m'  +  l  +  2  ^m 


g  ,  ,  a  .  b  ,  ,  =  0 
nn'  -1  m+1  m'  -1 


(19) 


Because  (15)  is  symmetric  with  respect  to  interchange  of  1  and 
2,  all  consequences  of  (19)  hold  as  well  when  a  and  b  are  inter¬ 
changed  . 

2 

Consider  possible  solutions  of  the  set  of  (2s+l)  equations 

(19)  . 
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m  =  m'  =  s  (^0 )  =$> 


a  b  =0 
s  s 


CASE  1  b  7^  0  =>  a  =  0 

-  s  s 


With  m'  =  Sj  (19)  reduces  to 


ms  a  b  +  •=•  g  g  ,  a  ,  b  ,  =  0 
ms  2  rs-l  m+1  s-1 


This  implies  =  0,  m  =  s  -  1,  s  -  2, 
bg_^  =  0.  With  m  =  1,  (19)  then  becomes 


g  g  ,  a  b  ,  ,  ,  =  0 
2  ro  rm'  o  m'  +  l 


(20) 

1  whether  or  not 


(21) 


and  with  m'  =  s  -  1  this  becomes  g  g  .a  b  =0  which  implies 

O  S  —  X  o  s 

a  =0.  Further  use  of  (20)  then  implies  a  =  0,  m  =  -1,-2,  .... 

-s.  We  therefore  conclude  for  CASE  1  that  there  is  no  solution 

with  b  ^  0,  a  =  0.  Naturally we  next  ask  if,  with  a  =  0, 
s  s  s 

there  is  a  solution  of  (19)  with  b  =  0  and  b  ,  ^  0.  If  not, 

we  then  ask  if  there  is  a  solution  with  b  =  b  ,  =  0,  b  -7^0. 

s  s-1  s-2 

And  so  on . 

CASE  2 

We  shall  proceed  by  treating  the  general  case:  we  show 
that  there  is  no  solution  of  (19)  with  b  =  b  ,  =  . . .  = 

■ -  1  —  b  b  J- 

b  >,1=0  a^d  b  ,  ^0,  k  =  1,2,  ...  ,2s .  To  prove  this  we 

S  •”  KH"  X  S  ~  KL 
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suppose  the  contrary  to  be  true.  With  m'  =  s-k,  (19)  becomes 


m  (s-k)  ambs_^.  +  2  ^m-l^s-k^-l^s-k+l  +  2  ^m^s-k-lanH-l^>s-k-l  ® 


(22) 


where  now  bg  =  0 . 


CASE  2.1  k  /  s,  k  /  2s.  The  argument  is  similar  to  that  of 

CASE  1.  Now  (22)  with  m  =  s  =>  a  =0,  then 

s 

”  m  =  s  -  1  =^>  as_2  =  then 

"  m  =  1  =>  3l  =  0 . 


Thus,  (22)  implies  am=  0 ,  m=  s,  s  -  1,  .  ..,  1.  With  m'  =  1 
in  (19)  we  again  obtain  (20)  which,  with  m'  +  1  =  s  -  k,  reduces 
to 


gogs-k-laobs-k  “  0 


(23) 


which  implies 


a  =  0 
o 


Then,  (22)  with  m  =  -1,-2,  .  .  . ,  -  s  ^  a  =  0,  m  =  -1,-2,  ...,-s. 

Thus,  all  the  a  =  0  and  NO  SOLUTION  is  possible  in  case  2.1 
m 

CASE  2.2  k  =  s 


(19)  with  m'  =  1  becomes  in  this  case 


g  g  a  ,  ,  b  =0 
2  ^m^o  mtl  o 


(24) 


9 


which  =>>  =  0 ,  m  =  s,  s  -  1,  .  .  .,  -s  +  1.  Then  (19)  with 

m  =  -s  +  1  =>  ... 

~2  ^-s^m'  a-s  ^m'  +  l  -  ® 

With  m'  =  -1  this  gives  g  g  ,  a  b  =0  which  a  =0. 

“  S  -1  "  s  o  —  s 

all  a  =  0  in  case  2.2  and  there  is  no  solution, 
m 


(25) 

Thus 


3 


CASE  2.3  k  =  2s 

(19)  with  m'  =  -s  +  1  in  this  case  ^ 


g  g  a  ..b  =0 

^m- s  m+1  -s 

With  ra  =  s  -  1,  s  -  2,  .  .  . ,  -s,  this  =^>  =  0 ,,  m  =  s,  s 

•  .  .  ,  —  s  +  1 . 

(19)  with  m'  =  -s  then  gives 


(26) 

1, 


-  ms  amb_s  =0  (27) 

with  m  =  -s  this  =>>  a  =0. 

-s 

Thus,  all  a  =  0  in  case  2.3.  and  there  is  no  solution.  This 
m 

exhausts  all  possibilities.  Thus,  we  conclude,  there  is  no  non¬ 
zero  solution  of  (19) .  In  other  words,  the  only  solution  of 
(15)  is  0^2  =  which  is  of  no  physical  interest. 

To  summarize  the  results  proved  so  far,  which  are  all  for 
N  =  2:  We  have  proved  that  Eq .  (3)  has  no  solution  with  ¥  ^  0 

when  X  =  0.  Since  from  Eq .  (4)  we  have  that  =  0  implies 
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X  =  0  if  Y  /  0,  the  only  solutions  with  Y  /  0  occur  for  ~o^  /  0. 
These  are  given  by  (14) . 

For  General  N: 

Take  the  scalar  product  of  (2)  with  0^ (3)0^ (4) . . .0^ (N)  and 
obtain 


’ S2  2  •  ( s*^+  2  ^j  '^k  ^1  ^^2  —  ^  0 ^  1 ) 0 2  ^ ^ ) 

k>  2  k>j>2 


(28) 


where  again  c..  =  <0 ^  |  s  1 0 ^>  . 

As  for  N  =  2,  we  consider  two  cases: 

2  "o’,  =  0.  Let  Z19  =  2  ct  •  If  2 

k>2  K  ^  k>2  K 


2  a,  ^  0  and  (b) 
k>2  K 

0,  the  eigenfunction  in 


(28)  is  well  known  to  be  Xg(l)Xs(2)  the  quantization  axis  being 

parallel  or  antiparallel  to  If  ^"12  =  (28)  reduces  to 

(3),  in  which  X  is  replaced  by  X  -  2  c.  -at  ,  so  that  again  the 

k>  j>2  3 

eigenfunction  of  (28)  is  X  (1)X  (2)  but  with  quantization  axis 

s  s 

undetermined . 

If  next  we  repeat  this  development  with  1  and  2  replaced  by 
i  and  j  we  obtain  instead  of  (28) 


+  2 
k>£ 

j 


ak 


0i(i)0j(j)  =  x  0i  (i)0  j  ( j ) 


(29) 
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with 


i: 


N 

=  2 

k=l 

Vi.  j 


o. 


(30) 


This,  of  course,  reduces  to  (28)  when  i  =  1  and  j  =  2.  The 

eigenfunctions,  determined  as  for  (28),  are  X  ( i ) X  (j)  where  the 

s  s 

quantization  axis  of  X.  is  arbitrary  if  2 • •  =  0  and  along  the 

s  13 

axis  of  2^j  if  2^j  0.  Now  let  i  =  1,  and  let  j  run  from  2  to 

N.  Then,  we  see  that  we  must  have  0 .  =  x„  for  all  j.  with  a 

3  s 

single  quantization  axis,  which  is  not  specified.  This  is  a 
consistent  solution,  since  2^j  =  2  =  (N-2)sz  for  all  pairs  i,j 

A 

with  z  the  quantization  axis,  a  consequence  of 

=  o’  =  <X£  Js1X£  >  =  sz,  independent  of  i.  Finally,  we  have 


X  = 


i<j 


-v  —  N(N-l)  2 

o  .  *o  .  =  — ^ s 

ID  2 


and  hence  from  (2),  a  =  Ns . 
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